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Abstract. Umbral calculus can be viewed as an abstract theory of 
the Heisenberg commutation relation [P, M] = 1. In ordinary quantum 

mechanics P is the derivative and M the coordinate operator. Here we 
shall realize P as a second order differential operator and M as a first 
order integral one. We show that this makes it possible to solve large 
classes of differential and integro-differential equations and to introduce 
new classes of orthogonal polynomials, related to Laguerre polynomials. 
These polynomials arc particularly well suited for describing so called 
flatenned beams in laser theory 



1. Introduction 

The purpose of this article is to construct a realization of the Heisenberg 
algebra in terms of a second order differential operator P and a first or- 
der integral one M. We then use this realization to construct families of 
orthogonal polynomials and study their properties. Finally we discuss ap- 
plications of these polynomials in mathematics (to solve integro-differential 
equations) , optics (to describe flattened beams in optics) and in astrophysics 
(to describe distortions of the microwave background radiation). 

This article is directly related to several research programs. The most 
general is that of umbral calculus, also known as finite operator calculus 
[1-4]. This has been described in many different ways. Here we view umbral 

calculus as an abstract theory of the Heisenberg relation P,M =1 and 
its mathematical implications. 

A related direction is that of monomiality, the essence of which is a sys- 
tematic investigation of the relations 

Pun= nUn-1, MUn = (1) 

leading to the representation of solutions of differential, difference and other 
equations as monomials Un = M'^uq [5-7] . 

The same algebra is of course obeyed by standard creation-annihilation 
operators in quantum physics. Monomiality and umbral calculus are closely 
related to the theory of Fock space ladder operators. For efficient recent 
applications of these concepts in mathematical physics see e. g. [5-11]. 

A recent article [12] was devoted to a systematic study of the realization 
of P as a difference operator and M as an operator the projection of which 
is the coordinate x. This was applied to show that continuous symmetries 
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like Lorentz or Galilei invariance can be implemented in quantum theories 
on lattices [12,13]. 

In this article we systematically investigate a different realization, namely 
one in which P is a second order linear differential operator in one real 
variable and M is a general first order linear integral operator. 

The problem is formulated mathematically in Section 2 where we also 
obtain the general form of the operators P, M in terms of one arbitrary 
function Y{x) and four constants. We obtain the boundary condition for 
functions u{x) defining the domain of the operators P and M. 

Section 3 is devoted to the basis functions Un = M'^uq, where uq is the 
seed or vacuum function. 

We show that we can either choose uq = 1 and then restrict the form of the 
function Y{x) or leave Y{x) arbitrary and choose uo{x) = Y'{x) {Y{x)+co)°' 
where cq and a are constants. We show that the basis functions Un{x) are 
eigenfunctions of a second order linear differential operator L with integer 
eigenvalues A„ = n + 1. The operator L is factorized in the sense that 
we have L = PM, where P, M are the differential and integral operators 
constructed in Section 2, respectively. 

In Section 4 we show that the eigenfunctions Un form an orthogonal set 
and obtain the corresponding measure and integration limits. 

The explicit form of the eigenfunctions is obtained in Section 5 for two 
different cases, namely uq = 1 and no(x) = Y'(x)Y(x)"', where a is a 
constant. In both cases the eigenfunctions are expressed in terms of Laguerre 
functions of non-trivial arguments. 

Section 6 is devoted to three types of applications: the solution of integro- 
differential equations, the description of flattened beams in optics and the 
distortion of background radiation in astrophysics. 

Conclusions and future outlook are presented in the final Section 7. 



(2) 



2. LINEAR DIFFERENTIAL AND INTEGRAL OPERATORS 
SATISFYING THE HEISENBERG ALGEBRA 

Let us consider two linear operators that we postulate to have the follow- 
ing form 

P = $2(x) Dl + + $o(x), 

M = g{x) i)-^ + k{x) 

where 

D = — 

^' (3) 
D~^u{x) = J u{x) dx ^ ^ 

Xo 

and $2, i,o(a^)) 9{x), k{x) are some sufficiently smooth functions, to be de- 
termined below. 

We impose the following requirements 

(1) The products P M and MP are differential operators 

(2) Their commutator satisfies the Heisenberg relation 



P,M 



1. 



(4) 
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The product P M will a priori contain a term proportional to D.^. ^ . The 
condition for it to be absent is 

$2(x) g"{x) + <^i{x) g'{x) + ^o{x) g{x) = 0, (5) 

where the primes denote derivatives with respect to the argument. 



The product MP and hence also the commutator 



P, M 



will apriori 



contain anti-derivatives ^of all orders. Indeed, the use of the Leibnitz 



formula extended to anti-derivatives yields 

oo 

b-\f{x)u{x)\ = Y,{-W{D^f{x))D-^^+^)u{x), 



(6) 



r=0 



where f{x) can be a function or an x-dependent operator. Calculating the 
product M Pwe obtain 

MPu = k{x) [<^>2{x) bl + ^i{x) + ^q{x)] u{x) + g{x) [$2(x) b,,- 
^'^{x) + ^i{x)] + Wi{x) - ^\{x) + ^o{x)]b-^ + ...\%^ 

(7) 

which will be a differential operators if 

^'i{x) - ^'^{x) + $o(x) = 0. (8) 

The coefficients of with k ^ 2 are all equal to (as differential conse- 
quences of eq. Q ) . 

The commutation relation Q then imposes two further equations 

2'^2{x)g'{x) + $'2(x)5(x) = 1 
k'{x) = 0, 

and a boundary condition at x = xq, namely 

[<^2{x)u{x) - <^'2{x)u{x) + ^>i(x)u(x)]U=^o = 0. 
Solving the first of eqs. (|9| either for ^2{x) or g{x), we obtain 



(9) 

(10) 



9{x) -- 
or 

^>2(x) 



+ 90 



J g{r])dr] + co 

Xo 



(11) 



go = 2 

where cq is an integration constant. Furthermore from the integration of the 
second of eqs. ([9| we find 

k{x) = ko (12) 
with ko a constant. Eqs. ^ and ([s]) provide expressions for ^i{x) and 
^o{x) in terms of $2(2^) and g{x). 

We shall use the second of eqs. ( pTj ) and get rid of the integral by intro- 
ducing a function Y(x), such that we have 



X 

g{x) = Y'{x), Y{x) = J giOd 



(13) 
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The (^-functions can then be expressed in terms of Y{x) and of some con- 
stants as 

1 



$2(2;) 
$o(x) 



y'(x)2 
1 

y'(x)3 
1 



[Yix)+co], 
[-3 {Y{x) + co)Y"{x) + (ci + 3) Y'ixf] 
+ co)iY"'{x)Y'{x) -3Y"ixf) 



(14) 



Y'{x)^ 
+{ci+3)Y'{xfY"ix)]. 



The boundary condition (10) can be rewritten in terms of the function Y{x) 
(see below). The operators P and M should be applied only to functions 
u{x) satisfying this boundary condition. 

The results of this section can be summed as a theorem 

Theorem 1. The operators 

P = Y^i^i^) + ^o]Dl + yj^s [(ci + 3)y'(x)2 - 3(y(x) + co)y"(x)] 

-FW [(^(^) + ^o){Y"'{x)Y'{x) -2,Y"{xf) + (ci + 3) Y'{xfY"{x)] , 
M = Y' b~^ + ko, 

satisfy the Heisenberg relation (M) and both P M and M P are second order 
differential operators. The function Y{x) satisfies the condition Y{xq) = 
and is otherwise arbitrary. The constants cq, ci, xq, are arbitrary. The 
operators P and M are defined for functions u{x) satisfying the boundary 
condition 

^ ;{(y(x)+co)y'(x)n'(x) + [(2 + ci)y'(x)2 (16) 



y'(x)3 

-(y(x) + co)y"(x)]n(x)} 



0. 



x=xo 



The inverse of the theorem is also true. All operators P, M satisfying the 
above properties are given by (15). 

3. MONOMIALITY AND BASIS FUNCTIONS 

3.1. General approach. The fundamental notion underlying monomiality 
is the existence of a sequence of basis functions ti„(x) satisfying 



Pun{x) = nu„_i(x), 
Mun{x) = n„+i(x), 

these two equations imply 

PMun{x) = (n + l)u„(x), 



P, M 



Un{x) = iun{x) 



The functions Un{x) are given explicitly as "monomials" 

Un{x) = M"uo(x), 



(17) 



(18) 



(19) 
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and in terms of the operator M and some "seed function" or " vacuum func- 
tion" tio(ic), which is to be defined. The use of the Heisenberg relation along 
with eq. ( 19 ) yields 

Pun{x) = nUn-i{x) + M''{Puo{x)). (20) 



To satisfy the first of eq. (17 1 for all n including n = 0, we must impose 

Pno(x) = 0. (21) 

The above conditions can be viewed in at least two different ways 

a: As a condition on the operator P and thus on the function Y{x) 
b: As a condition on the seed function uo{x). 

It is evident that uq{x) plays the role of a physical vacuum. 

3.2. Conditions on the operatorP for a constant seed function. Let 

us put uq = 1 (or any other non-zero constant). Condition (21) implies 
^o{x) = that is, in view of eq. (15) 

{Y{x) + Co) {Y'{x) Y"'{x) - 3 Y"{xf) + (ci + 2,)Y' {xfY" {x) = 0. (22) 

The function Y{x) is therefore no longer arbitrary, but it is subject to eq. 
(22). This equation has a three-dimensional Lie point symmetry group, 
generated by the Lie algebra 

Xi = dx, 

X2 = xd„ (23) 
X3 = iY + Co) dr 

which can be used to reduce eq. ( p2| ) to quadratures. However we obtain 
implicit solutions of little use in the present context. We therefore use an 
alternative approach for uq = 1, namely we return to the equations solved 
in Section 2 and impose ^o{x) = from the beginning. Eqs. (|5] [s]) with 
^o{x) = imply 



$i(x) = ^>2(x) + a, 

1-2/3 
<I>^(z)-2a' 



where a, (3 are integration constants. Eq. ^ implies that <I>2(2;)satisfies the 
condition 

(1-2/3) $2(x) $2(x) + /3 [$'2(x)] ^ - a (1 + 2 /?) ^'^{x) + 2 = 0. (25) 

By solving the above equation for a = we obtain a simple but interesting 
solution 

$2(x) = ^^±g^, cl>,(^)=cl>'^(^) = _l_zi^ (^ + ^)-'^, <I>0(x)=0, 



1+9- 



(26) 



Eqs. (14) then imply that 



Y'{x) = A{x + ^y, co = ^(xo+7)''+\ 

n^) = <FT[(^ + 7)'^'-(^o + 7)^+^], ci = -2±f. 
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The operator P reduces to 



P 



A(q + 1 



(28) 



and is self-adjoint. 



The boundary condition (16) is satisfied identicaUy for u[x) = const once 



the constants cq i are those given in eq. (27). 



3.3. Condition on the seed function u{x) for arbitrary Y{x). Let us 

now consider the operator P of eq. (15) with Y{x) (and hence g{x) = Y'{x)) 
arbitrary. The seed function must satisfy eq. ( |21[ ). From eq. ([5]) we see that 
uq{x) = Y'{x) is a solution of eq. (21). A second linearly independent 



solution is easily obtained using the Wronskian. The general solution of eq. 



(21) can be written as 

Y'{x) [ai + a2{Y{x) + cq) 



Uq{x) 



-ci-21 



Ci / -2, 



(29) 



or 



uo{x) = Y'{x)[ai + 02 ln(y(x) + co)] , ci = -2. (30) 

Since we are dealing with linear equations we can consider separately the 
cases oi = 1, 02 = and oi = 0, 02 = 1 when imposing the boundary 
condition (16). The boundary condition is satisfied for (2;) = Y' {x)(Y{x) + 
co)"'^^"^ for all values of ci, in particular for the case ci = —2, i.e. uq{x) = 
Y'{x). It is not satisfied for the term with a logarithm, so we discard the 



solution (30). 



3.4. The eigenvalue problem for the basis functions Un{x). The func- 
tions Un{x) defined by the relations (17 18 19) with uo{x) satisfying eq. 
(21 ), will be eigenfunctions of the linear operator L = P M and the situation 



can be summed up according to the following theorem 

Theorem 2. The functions Un{x) defined by 

Un{x) = M"-uo{x), 
Puo{x)=0, 

satisfy the second order ODE 



LUn{x) = (n + l) Un{x), 

L = /2(x) Dl + f^{x) D, + fo{x) 



(31) 



(32) 



with 

f2{x) 

fo{x) 



y'{x)3 



yrj^[Y{x) + Co], 

{ [Y{x) + CO + Uci + 3)] y'(x)2 -SkoiYix) + Co) Y"{x)}, 

ko{Y{x) + Co) {Y'{x) Y"'{x) -3Y"{xf)+, ^ 
+[A:o(ci + 3) + Y{x) + cq] Y'{xfY"{x)- 
-(ci + 3)r(x)4 

(33) 



'y'(x)4 



as long as Y(x) satisfies the boundary condition (16). The operators M , P 
are those of Theorem^ The seed function uq{x) must be chosen as in eq. 
( 29) for Y{x) arbitrary, or Uo{x) = 1 for Y{x) as in eq. (21). 
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The seed function uo{x) satisfies eq. (16), but this does not guarantee the 
same for all Ui{x), i ^ 1. 

We will see in Section 5 that this imposes a condition on the constants 
co,i- 



3.5. Eigenfunctions as functions of two variables. We have constructed 
the monomials m„(x) = M"mo as functions of one variable x. They also de- 
pend in a significant way on the parameter ko. Let us, for the purposes of 
this paragraph, change notation and write 



Trnix,y) = Unix), y = ko. 



(34) 



We then have 



M = Yix)D-^ + y. 



(35) 



The functions 7r„(x,y), as functions of two variables, satisfy a partial 
differential equation 



Indeed, we have 



dy 

d7Tn{x,y) 
dy 



pTTnix,y). 



(36) 



dy 



We have dM/dy = 1, and eq. (17) leads to eq. (36). Eq. (36) can be 
written explicitly as a second order partial differential equation, using eq. 
(15). In particular, for a constant seed function no = 1 we have P as in eq. 
(28) and eq. (36) reduces to 

dTTn{x,y) ] 



dy 



A{1 + . 



-D^X^ ''D^TTn, q / -1, 



(37) 



where we have put 7 = 0. 

Eq. (37) is a linear heat equation with variable conductivity x^^'^ and the 
monomials n„(x) = 7rn{x,y) = M^uq are solutions of this equation. 

Let us put 



y = A{l + q)t, 1 



(iV/2). 



Eq.(37) now is 



du{x, t) 
dt 



DxX^ Dxu{x, t). 



(38) 



Its Lie point symmetry algebra can be calculated using standard methods 
[14]. For AT = and = I the algebra is six dimensional and the N = ^ 
case is isomorphic to the = one, i.e. to the symmetry algebra of the 
constant coefficient heat equation. A basis for the symmetry algebra for 
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= and iV 

Xi = 

X2 = 

X3 = 

X4 = 

X5 = 

Xe = 



3 IS 

dt, 



(39) 



1 



_ ^ 1 
N 2(2 - N) 



udu, 



N 



txdr 



1 



„2-N 



(2-iV)2' 



+ 



N 



t] udu, 



Ox - —X 2 udu, 



JV 

tX 2 dr 



rN iV-2 

[-tX 2 + 



:X 2 ]udu. 



N 



udu- 



The fact that the two algebras are isomorphic suggests that the two equa- 
tions could be transformed into each other by a point transformation (it is 
a necessary condition for the existence of such a transformation). This is 
indeed the case here. We put 



y = t, z = 3x^^^, w{z,y) = x^^^u{x,t). 
Then, if u[x, t) satisfies 



du d 
dt dx^ 



,4/3 



du 
dx^ 



w{z,y) will satisfy 



and vice versa. 



(40) 
(41) 
(42) 



For N ^ 0, ^,2 the symmetry algebra of eq, (38) is four-dimensional with 
basis given by Xi, X2, X3 and Xq of eq. ( 39 ) . The symmetry group in this 
general case is GL{2,TZ) and the equation (38) can not be transformed into 
the usual heat equation ( 42 ) . 



4. ORTHOGONALITY PROPERTIES OF THE 
EIGENFUNCTIONS 

Let us first recall some well known results from the theory of linear oper- 
ators [?, 14] The adjoint L"*" of the second linear differential operator L of 
eq. ( 32 ) is defined by the relation 

b b 

J U2{x) L Ui{x) dx = j Uiix) L+U2{x) dx, (43) 

a a 

where U2,i{x) satisfy the boundary condition, 

{U2{x) f2{x) U[{x) - Ui{x) (/2(x) Ulix))' + /i(x) Ui{x) U2{x)) t = 0, 

(44) 

and we have 

L+ = f2{x) bl + (2 f'^[x) - h[x)) bx + f^- f[{x) + fo{x). (45) 
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Let US introduce the eigenfunctions of L and 



LUn{x) = XnUn{x), 
L+Vm{x) = \mVm{x). 



(46) 



The functions Un{x), Vm{x) form mutuahy orthogonal sets 
b 

Vm{x) Un{x) dx = 0, m u Am / A„, (47) 



where Vmix) (in the domain of L^) and Un{x) (in the domain of L) must 
satisfy the same boundary conditions ( |44[ ) (for some a and b ). Now let us 
request that the eigenfunctions Vn{x) be proportional to Un{x) 

Vn{x) = w{x) Un{x). (48) 

This implies that w{x) must satisfy the following first order ODE 

f2{x) w'{x) + (/^(x) - /i(x)) w{x) = 0. (49) 

The boundary condition for Ui{x) = n„(x)and U2{x) = Vm{x) = w{x) Um{x) 
reduces to 

[f2{x) w{x){Um{x) <(x) - U^{x) Un{x))] ^ = 0. (50) 



If w{x) satisfies ( [49| ), ( [50| ) then the eigenfunctions Un{x) of the original 
operator L will be orthogonal with the weight w{x) 

b 

w{x)Umix)Un{x) dx = 0, Xm ^ ^n- (51) 



Applying the above general results to the operator L of ([2| we obtain the 
following theorem 



Theorem 3. The eigenfunctions Un{x) of eq. (32) will satisfy the orthogo- 
nality relation 

h 



^ w{x) Urn{x) Un{x) dx = NnmSnm, (52) 
a 

with 

{Y{x) + Co) '^+^ ^Y{x) 

w^^) = ^F7(^^ 6^=0 \ (53) 



provided they satisfy the boundary conditions 

^^^^^ (iim(xX(a;) - n^(a;)n„,(x)) I =0. (54) 



(y(x)+Co)^i+3^^y( 



\Y'{x)\ 



We mention that the boundary condition (54) for orthogonality and (16) 



for ([1| [2j) to hold must be both satisfied and that xq, a, b are apriori inde- 
pendent constants (to be determined below). 
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5. EXPLICIT FORMS OF THE EIGENFUNCTIONS 

5.1. Eigenfunctions for Uo{x) = Y'{x) [Y{x) + cq] " with y(x)arbitrary. 
Here we proceed in tlie spirit of Section (3.3) and put 

ci = -2-a, -uoCx) = y'(x) + Co) " (55) 



i.e. we start from the second solution in eq. ( 29 ) . The first one in ( 29 ) is just 
the special case of a = 0. The condition (21 ) is satisfied as is the boundary 

ext e 

"(y(x) + co)"+i 



condition (16) for uq(x). The next eigenfunction calculated according to eq. 
((191) is 



ui{x) = Y'{x) 



a + 1 



+ ko{Y{x) + Co) 



1 

a + 1 



Substituting into the boundary condition (16) we obtain Cq"*"^ = 0. It follows 
that we have 

Co = 0, a > -1. (57) 

Using eq. (19) we can easily generate U2{x), U3{x), .... Inspired by the form 
of these functions we make the Ansatz 



(56) 



Unix) = Y'{x) y(x)"5n(- 



y(x). 



(58) 



and substitute it into the eigenvalue equation ( |32[ ) . We find that the function 
gn satisfies the generalised, or associated [16, 17] Laguerre equation 



. Y{x) „ 

( 7— ffn + 

fcn 



a + 1 - (- 



y(x). 



gjx) +ngn{x) 







(59) 



Thus eq. (32 ) is solved in terms of generalised Laguerre polynomials lI^\z). 
Substituting (58) into the boundary condition (16), we see that (16) is sat- 
isfied for all values of n. The boundary condition (54) for orthogonality is 
satisfied for all n and m if we choose 



= 0, b = oo, 



Y{x) 



-oo 



(60) 



(we have put, with no loss of generality, xq = 0). Eq. (31) gives us the 
explicit form of the eigenfunctions for all n. 

Let us sum up the previous results as the following theorem 

Theorem 4. The eigenvalue problem 

y^y(x)<(x) + yrf^ {[Y{x) + ko{l - a)] Y'{x)^ - 3koY {x)Y" (x)} u^x) 

-yT^jWx) {Y'{x)Y"'{x)-:iY"{xf) + [ko{l-a)+Y{x)] Y'{xfY"{x) 

-(1 - a) y'(x)4| Un = {n + 1) Unix), 

(61) 

with the boundary condition 
1 



(62) 
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is solved by the monomials 

Un{x) = M'^uoix), uo{x) = Y'{x)Y{x)'' , M = Y'{x)D-^ + ko. (63) 
Explicitly the solutions are 



0=0 



Unix) = Y'{x)Y{xrY.]^Xn-" )Y{xyK~' 



(64) 



y'(x)y(x)"L(")(-^), a 



> -1, 



where L^"^ {z)is a generalized Laguerre polynomial. The function Y{x) and 
the constants k^, a satisfy 

Y[x) 



lim 



-oo, a > — 1. 



The orthogonality relation is 



Y{x) " ^Y{x) . s . s , f 1 \ax 
— — — 6*^0 'Un{x) Um{x) dx = {-korSr, 

\Y'{x)\ 



(65) 



(66) 



5.2. Eigenfunctions for uq(x) = 1, Y'(x) = A {x + 'j) Let us now take 
Y{x), Co and ci as in eq. (27). The eigenvalue problem p2|) reduces to 



(x + j)'-'^Dlunix) + ^^{A (x + 7) - Hq - l)(x + j)-i}D,Un{x) 
= nun{x), 

_ (67) 

(68) 



the boundary condition (16) reduces to 



(a; + 7)^~V(x)U=^„ =0. 



Condition (68) is satisfied for uq{x) = 1. 
We have 



ui{x) = Muo{x) = A{x + jy{x-xo) + ko 



and eq. (68) for tii(x)implies 



7 = -xo, Co = 0, 



(69) 
(70) 



so that we find 
P ' 



A(q+1 



a,{x-xoy-'iD,, M = A[x-XQfb-^ + kQ. (71) 



We can easily implement eq. ( 19 ) and get 



w^E(;)*r; 



A\x - xo)^(«+^) 



{q + 2){2q + 2,)...{{j-l)q + j) 



(72) 



Substituting u„ into the boundary condition ( |68[ ) with 7 = —xq we obtain 
the condition 



q>-l. 



(73) 
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Eq (72) allows us to relate u„(x) to the generalised Laguerre polynomials 



according to the identity 

Un{x) 



a 



Z = 



^0 ^ (^) 



(74) 



where we have set xq 



ko (q+l) • 

0, with no loss of generality. 



The substitution (74) reduces the eigenvalue problem (67) for 7 = —xq = 
to the generalized Laguerre equation for L^\z). The orthogonality rela- 



tion ( 52 ) in this case specializes to 



A 
q + l 



n 



9 

9+1 



1 - 



1 Ax1+'^ 

-e''0(i+^'> Unix) Umix) 



(75) 



A 



2q_ 

9 + 1 1 



(-A;o)'+i(5n, 



In conclusion we state the following theorem. 
Theorem 5. The eigenvalue problem 



<{^)+ i^x'i + {l-q)x'^ <(x) 



A(q+l) 
ko 



n x'^ ^Un{x) = 0, 



u'nix)\ 



x=0, 



> -1, 



(76) 



is solved by the monomials 



(77) 



Unix) = M"l, 

M = AxW-^ + ko. 

The eigenf unctions Un{x) are expressed in terms of generalised Laguerre 
polynomials by eq. (74)- Their normalization is given in eq. 



75). 



As a further remark let us note that the polynomials we have derived 
cannot be framed within the context of the Sheffer type family, the operator 
P is indeed a function either of x and of the ordinary derivative (see 
ref. [7] for further comments). More generally they belong to the Laguerre 
family which provides polynomial forms whose generating function cannot 
be expressed in terms of exponential functions. 

6. APPLICATIONS 

6.1. Applications to the solution of integro-difFerential equations. 

An important aspect of umbral calculus is the umbral correspondence [1-4]: 
the correspondence between results obtained in different realizations of the 
Heisenberg algebra. 

Let us first consider a simple example, namely the initial value problem 
for a first order partial differential equation 



dy_ 
dr 



d 

(^+a:)y, y{0,x) 



I. 



(78) 



It can be solved by the method of characteristics, or in a more formal way 
by putting 



y{T,x) 



(79) 
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Using the Baker-Campbell-Hausdorff formula [18], which in this case re- 
duces to 



we obtain 



y{T,x) 



^ kl 



(80) 



(81) 



fc=0 



Now let us replace the PDE ( 78 ) by an operator equation 

yr{T,M) = {P + M)y{T,M), (82) 

where P and M are the operators ^ studied in the previous sections. In eq. 
(82 ) y(T, M) is an operator function. We will apply both sides of eq. (82 ) to 



the seed function uq and this will turn eq. ( 82 ) into an integro-differential 
equation for a function f{T,x). 



A formal solution of the operator equation (82) is obtained from the 



expansion (81) simply by inserting the operator M instead of x. We ap- 



ply this operator y{T,M) to the seed function uo{x) and obtain f{T,x) 
y{T,M)uo{x). i.e. 



/(r,x) 



oo ^ 



(83) 



A:=0 



where Uk{x) are the basis functions of Section 5, ultimately expressed in 
terms of generalized Laguerre polynomials. 

The function /(r, x) is at least a formal solution of the equation 

yr{T,M)uo{x) = {P + M)y{T,M)uo{x), (84) 

in our case the integro-differential equation 

fr{T,x) = <!>2Dlf{T,x) + <^M{T,x) + {<!>o + ko)f{T,x) (85) 

/•X 

+Y / f{T,x)dx. 
Jo 

Here ^2, ^1 and <I>o and Y are the functions determined in Section 2. The 



formal solution (83) is a real solution of eq. (85) if the series converges. 



Let us sum up this example. We start from the linear partial diffeential 



equations (78) for which we know the solution (81) of a given initial value 



problem. We associate an operator equation (82) and an integro-differential 



equation (85 ) to (78) via the umbral correspondence. A solution of an initial 



value problem for (85 ) is obtained from (81 ) by the umbral correspondence: 



we replace the powers x" by the basis functions Un{x) for the appropriate 
operators P and M of Sections 3, 4 and 5. The solution (83) corresponds 
to the initial value 



/(0,x) = uo{x) 



(86) 



for eq. ( 85 ) . 



More generally, let us consider an evolution equation of the form 

yr{T,x) = F{d.,,x)y{T,x), y(0,x) = /(x), (87) 
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where F(dx,x) is a polynomial in dx with coefficients that are power series 
(or polynomials) in x. A solution of the initial value problem (87) is given 
by 



y{r,x) 



where (88 ) can be written (at least in principle) as a power series by applying 
the Baker~Campbell~HausdorfF formulas in eq. (88): 

y{T, x) 



j=0 



T]X-' 



(or the power series ( 89 ) can be obtained in any other way) . 
To the function (89) we associate another function, namely 



/(r,x) 



y(r, M)uo = ^ Cj(r)'Uj(x). 

j=0 



(90) 



The function ( 90 ) will be a solution of the Integra-differential equation 

yr{T, M)uo = F{P, M)y{T, M)uo, (91) 

i.e 

fr{T,x)=F{P,M)f{T,x), (92) 



with initial condition 



f{0,x) = Y,Cj{0)u,{x). 

3=0 



(93) 



In other words, if we know how to solve the PDE (87), we know how to 
solve the integro-differential equation ( 92 ) . 

6.2. Description of flattened optical beams. In this subsection we will 
emphasize the extreme usefulness of the above family of orthogonal polyno- 
mials in applications to optics. 

We define the orthogonal function 



<^n{x,y) = ane ^m.+d"" Vi'?)(x, y), 
with Q„ a suitable normalization constant such that 



^n{x, y) (^m{x, y) dx = 6n 



(94) 



(95) 



The functions TTn\x,y) are defined similarly as in Section 3.5, i.e. we put 

7ri^\x,y) = {Axm-'+yri, y = ko- (96) 

An idea of the shape of the above family of orthogonal functions is given in 
Fig. 1 where we have plotted the ffi'st three "modes" for g = 3. The shape 
is that of the so called flattened optical beams. These beams are different 
from the usual Gaussian beams since they have flat transverse distribution 
as shown in Fig. 2 and Fig. 3. The advantage offered by flattened beams 
with respect to ordinary Gaussians is that they provide larger suppression of 
thermal noise on the mirror surface because of a better average on the surface 
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Figure 1. The first three flattened q = 3 beam modes nor- 
mahsed to unity (n = continuous, n = 1 dot, n = 2 dash). 
We take A = 1 and y = —1. 



Gaussian B«am FlaLlop B^am 




Figure 2. Transverse distribution of a Gaussian and a Flat- 
tened beam, the difference in intensities ( denotes the beam 
waist, with L being the length of the cavity and k the wave 
vector of the optical fleld ) 




Figure 3. Transverse distribution of a higher order flat- 
tened mode 

fluctuations, they are therefore particularly useful for high power lasers. It 
is evident that if we take a transverse section of the flat top distribution in 
Figs. 2, 3 we get the distribution similar to that given in Fig. 1 for n = 
and n = 1 respectively [19] used in optics to treat a light beam whose cross 
section has an intensity as uniform as possible. 
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■11 




t-, 



Figure 4. Flattened beam mode distribution q = 3 (nor- 
malized to unity) and comparison with the fundamental, a) 
n = 4 ,b) n = 30 



A typical example of flattened beam is a super-gaussian 

fix) oc e--" (97) 

whose profile becomes more and more flat as p increases. 

It is evident that the function <^„(x,y) is a super-gaussian for n = and 
that for n > we have higher order super-gaussian flattened modes. 

The distribution of a higher order modes can be obtained from the square 



moduli of the functions (94) and are shown in Fig. 4 for the cases n = 4 
and n = 30. 

The advantage offered by this family of orthogonal functions is two-fold 

a: They provide the natural set for the expansion of flattened beams 
b: They offer the possibility of treating the higher order modes and 
not only the fundamental one. 

The study of the propagation of the above family of flattened beams can be 
performed using the set of functions introduced in this article. We postpone 
this to a forthcoming investigation. 
It is worth pointing out that we have 

c^„(x,y) = a„(-y)"e-i4")(z) (98) 

and we can hence discuss the relevant evolution using z as the independent 
transverse coordinate. 

It is also worth stressing that the operators P, M can be used to form 
other Lie algebras, different from the Heisenberg one, similarly as x and 
p = —i/dx are used to form e.g. su(l,l). In turn these are well suited to 
describe the optical cavity elements and filters exploited to flatten the beam 
distribution, as it will be shown in a forthcoming investigation. 

6.3. Applications in Astrophysics. The Sunyaev-Zeldovich effect [21- 
23] is the distorsion of the cosmic microwave background radiation spectrum 
by the inverse Compton scattering of high energy electrons. When describing 
this effect the authors [21,22] obtained the diffusion equation (the Sunyaev- 
Zeldovich equation): 
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Eq. (99) has a symmetry algebra with basis: 



^5 

^6 



rdr + 1^ log - ^ [9r + 3 log e + l] wd,,, 

T^dr + log - ^ [9r2 + 2r + 6r log ^ + log^ (\ wd^, 



(100) 



-w 



[3r + log (\ wdu 



This Lie algebra is isomorphic to the algebra given in eq. ( 39 ) , implying 



that the Sunyaev-Zeldovich equation (99) might be equivalent to the heat 



equation. Indeed it is and the equivalence is realized by the transformation 

t = r, x = log^, n(x,t) = C^/2e3^uj(C,T), (101) 

(see also [24]) If w{^^t) satisfies eq. ( [99] ) then u{x,t) satisfies the heat 
equation ut = Uxx- Hence the Sunyaev-Zeldovich equation can be solved 
in terms of the functions (72) constructed in this article (for q = 1 and 
ko = 2 At). 

7. CONCLUSIONS 

One way of summing up the results of the present article is the following. 
We have constructed the most general operators P, M of the form ^ that 
satisfy the Heisenberg relation (|4]) and used them to construct the monomials 
Un{x) = M"tio(x). They are eigenfunctions of the linear operators L = PM 
corresponding to positive integer eigenvalues Xn = n + 1 . 

The second order operator L given in eq. (32) is thus factored into a 



product of 2 operators P and M. The usual factorization of a differential 
operator is into two lower order differential operators [25,26]. Ours is highly 
non standard: a differential operator P times an integral one M. 

The functions Un{x) are expressed in terms of generalized Laguerre poly- 
nomials lIt'\z), or equivalently , the confluent hypergeometric function 



2' 



or a 



2Fo{a,b; z). This includes the Hermite polynomials for a 
but none of the other classical orthogonal polynomials, related to the hy- 
pergeometric function (rather then the confluent one). 

There is a good reason for this. We have imposed the form ^ on the 
operators P, M and the monomiality condition then leads to an eigenvalue 
problem in which the order of the polynomials n enters in the eigenvalues 
only and enters linearly. For all other classical orthogonal polynomials the 
dependence on n is more general. 

We are looking into possible generalizations in order to obtain monomial 
realizations of other classes of orthogonal polynomials. This can e. g. be 
done by imposing all the properties of monomiality, but allowing P, M to 
depend on n, 

PnUn{x) = nn„_i(x), MnUn{x) = U„+i(x). (102) 
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The Heisenberg relation should in this case be modified to 

{Pn+lMn - Mn-lPn) = Un[x) (103) 

As also stressed in the previous section the possibility of introducing an 
extra variable makes it possible to consider our polynomials (or functions) 
as depending on two variables (and also on two parameters n, a). 

A further topic that is being pursued is that of applications, as outlined 



in Section 6. A very special case of eq. (92) was studied in Ref. [8]. We 



intend to study eq. (92) in general with P, M as specified in this article, for 
a wide choice of the functions F{P^ M). 

Regarding applications we have stressed that the form of the introduced 
orthogonal polynomials is ideally suited for the study of optical flattened 
beams. 
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